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We study the twist-4 contribution to the unpolarized nucleon structure functions, 
Fl and F2, in an approach where the vacuum structure of QCD is described 
by a dilute medium of instantons (p/R <C 1). The dominant power corrections 
come from the twist-4 quark-gluon operators and are of the order l/(p 2 Q 2 ). The 
contributions from four— quark operators are suppressed. The approach predicts 
large l/Q 2 -corrections to Fl, in agreement with the results of QCD fits to the 
data. 



The operator product expansion of QCD explains not only the scaling 
behavior of structure functions at asymptotically large Q 2 — including the 
logarithmic corrections, which can be computed using perturbation theory — 
but allows to classify also the power (1/Q 2 -) corrections to scaling. The latter 
are determined by the expectation values of certain operators of non-leading 
twist (> 2) in the target hadron state Ll While a fully quantitative theory 
combining logarithmic and power corrections is presently out of reach, it is 
interesting to study these higher-twist matrix elements in their own right. 
These are "universal" measures of the non-perturbative correlations of the 
quark and gluon fields inside hadrons, which appear also in other contexts, 
e.g. in the-description of "intrinsic" heavy quark contributions to light hadron 
properties El. Knowledge of these matrix elements would be useful for a more 
precise extraction of the leading-twist parton distributions from the DIS data, 
and for a better understanding of parton-hadron duality!! 

Two basic types of twist-4 operators appear in the l/Q 2 -power corrections 
to tlio-second moments of the unpolarized structure functions (for details, see 
RcfsEJu). One are quark-gluon operators, which describe the influence in DIS 
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of the non-perturbative gluon field in the target on the "struck" quark: 

(p\tp f j a (D^F^)ij f \p) - trace = 2A f 

(p\ $n J 75i {v a F Pl + F ai V^ ipf \p) - trace = 2 5/ 

(1) 

Here, F a p [F a p) denote the (dual) gluon field. The other are four-quark oper- 
ators, describing the interference contribution to DIS from scattering off two 
different quarks in the target: 

9 2 (Pi ipf y7a75^/ ipf y 7/375^/' \p) = ZCffi [p a Pf3 - ~^9aA ■ (2) 

In this context two questions arise: 

• What are the relevant non-perturbative scales determining the values of 
the twist-4 matrix elements, A, B and CI 

• Are the matrix elements of both types of operators of equal magnitude 
or is one of them dominant? 

It is interesting to address these questions from the point of view of the 
instanton vacuum, where one studies non-perturbative effects generated by a 
dilute "medium" of instantons and antiinstantons. These vacuum fluctuations 
are known to be responsible for the spontaneous breaking of the chiral sym- 
metry of QCD and a variety of other effects; seeu for a review. The properties 
of the instanton "medium!! have been derived within a variational approach 
to the QCELground states, and have also been extensively studied in lattice 
simulations a. A crucial property is that the average size of the instantons in 
the medium, p = (600 MeV) -1 , is small compared to their average distance, R, 
with p/R ~ 1/3 ( "dilutencss" ) . This ratio provides a small parameter which 
allows to classify non-perturbative effects generated by the instanton medium. 
For instance, the dynamical quark mass, M, which appears in the spontaneous 
breaking of chiral symmetry, is parametrically of the order M ~ (p/R) 2 p~ x . 
As a result, e.g. the pion weak decay constant is parametrically suppressed: 

fc 2 <p- 2 

(3) 

Similar parametric estimates can be performed for other observables. 

Following the same philosophy we can estimate also the magnitude of the 
twist-4 nucleon matrix elements, (01) and (0)0'cl When computing the matrix 
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element of the twist-4 quark-gluon operators, ([!]), the gauge field in the oper- 
ator is replaced by the field of a single (anti— ) instanton, which then interacts 
with the fermion fields in the operator through its fermionic zero modes. (The 
same zero-mode induced interaction is instrumental in the dynamical breaking 
of chiral symmetry.) One finds that the isoscalar nucleon matrix element is 
determined by an integral of the form 

fc 2 <p- 2 

Contrary to the pion decay constant, (^), this integral is quadratically di- 
vergent, so that its value is determined by the square of the ultraviolet cut- 
off — the inverse instanton size, p^ 1 . A more quantitative estimate gives 
B u + Bd ~ 0.9 JT 2 — (570 MeV) 2 . The matrix element has a large numerical 
value, because of the smallness of the average instanton size. The correspond- 
ing matrix element of the first operator, (Q), is accidentally suppressed, since 
D 1 F 1 p = in the field of one instanton. For the matrix elements of the 
four-fermionic operators, on the other hand, we find 

C uu ,C dd ,C u d ~ F 2 ~ (|) * log (£) p-\ (5) 

These matrix elements are parametrically suppressed relative to that of the 
operator B, ([|). In addition, they are suppressed in the coupling constant, g 2 . 
Numerically, a value of (50 MeV) 2 for the isoscalar matrix element was found 
in RefQ — two orders of magnitude smaller than for the quark-gluon operator! 

To summarize, the parametric order of the higher-twist matrix elements 
obtained from the instanton vacuum implies that: 

• Matrix elements of twist-4 four-quark operators, (|^), are parametrically 
(and numerically) small compared to those of quark-gluon operators, ([!]). 

• The relevant scale for twist-4 quark-gluon operators is the square of 
the average instanton size, p~ 2 = (600 MeV) 2 . The dominant power 
corrections are thus expected to be of the order l/(p 2 Q 2 ). 

These findings are in agreement with the results of QCD fits to the data fox 
Fl and F2 , which show a sizable 1/Q 2 corrections to Fl consistent with (||) . Q 
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